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Abstract We study charmonium (ccnfi), bottomonium (bbnn) and exotic (ccnn and bbnn) four-quark states 
by means of a standard non-relativistic quark potential model. We look for possible bound states. Among 
them we are able to distinguish between meson-meson molecules and compact four-quark states. 
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1 Introduction 

The discoveries on several fronts of unusual char- 
monium states like A(3872) and Y(4260) and open- 
charm mesons with unexpected masses like Z?* 7 (2317) 
and Z?2(2308), have re- invigorated the study of the 
hadron spectra. Their anomalous nature has trig- 
gered several interpretations, among them, the exis- 
tence of compact four-quark states or meson-meson 
molecules. This challenging situation resembles the 
long-standing problem of the light-scalar mesons, 
where it has been suggested that some resonances 
may not be ordinary qq states, though there is lit- 
tle agreement on what they actually are ^. In this 
case, four-quark states have been justified to coexist 
with qq states because they can couple to J PC = ++ 
without orbital excitation ® . 

Four-quark systems present a richer color struc- 
ture than standard baryons or mesons. Although the 
color wave function for ordinary mesons and baryons 
leads to a single vector, working with four-quark 
states there are different vectors driving to a singlet 
color state out of colorless or colored quark-antiquark 
two-body components. Thus, in dealing with four- 



quark states an important question is whether we 
are in front of a colorless meson-meson molecule or 
a compact state (i.e., a system with two-body col- 
ored components). Whereas the first structure would 
be natural in the naive quark model, the second one 
would open a new area in hadron spectroscopy. 

We have derived the necessary formalism to eval- 
uate the physical channels probability(singlet-singlet 
color states) in an arbitrary four-quark wave func- 
tion. For this purpose we have expanded any hidden- 
color vector of the four-quark state color basis, i.e., 
vectors with non-singlet internal color couplings, in 
terms of singlet-singlet color vectors. Such a proce- 
dure gives rise to a wave function expanded in terms 
of color singlet-singlet nonorthogonal vectors, where 
the determination of the probability of physical chan- 
nels becomes cumbersome. However, it allows to dif- 
ferentiate among unbound, compact and molecular 
four-quark states. 

2 Formalism 

There are three different ways of coupling two 
quarks and two antiquarks into a colorless state: 



[(«i«b)(«3?4)] = {|3i23 3 4),|6 12 6 3 4)} = {|33)f,|66)f} 

[(ft&Xfcfc)] = {|ll3l24>,|8 13 8 24 }} = {|ll) c ,|88) c } 

[(gi&Xfcfs)] = {|li4l 2 3),|8 14 8 23 )} = {|lT') c ,|8'8') c }, 



(la) 
(lb) 
(lc) 
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each forms an ortho-normal basis. Each coupling 
scheme allows us to define a color basis where the 
four-quark problem can be solved. The first basis, 
Eq. (Taj) , being the most suitable one to deal with the 
Pauli principle is made entirely of vectors containing 
hidden-color components. The other two, Eqs. (|lb[) 
and (fTcjl . are hybrid bases containing singlet -singlet 
(physical) and octet-octet (hidden-color) vectors. 

An arbitrary four-quark wave function can be 
expressed in terms of physical components, singlet- 
singlet color states, by means of an infinite expansion 
than can be resumed as ^ , 



l*> = \{pQ+Qp) i 1 2 I*) 

2 V / 1 — cos 2 a 

+ 1(pQ+Qp)i — ^1*). (2) 

2 V / 1 — cos 2 a 



Thus, one obtains two hermitian operators that are 
well-defined projectors on the two physical singlet- 



singlet color states 

V\u }c = (r<) <)p 



1 



2(1 — cos 2 a) 

= ( pQ+QP >2(do^)- < 3 > 



Thus, given an arbitrary state its projection on 
a particular subspace E is given by \^)\ E = P E \^>). 
Then, the probability of finding such an state on this 
subspace is 

B |<*|*)| B - (*|PlP B |*> = <*|i^|*) = {9\P B \9) . 

(4) 

Therefore, once the projection operators have been 
constructed [Eq. ©], the probabilities for finding 
singlet-singlet components are given by, 



P"*>([11]) = <*|7>|u>J*> 
P'*>([1'1']) = (¥|7Vi'>.l*) 



(5) 



Using Eq. © it can be easily checked that pl*^ ([11])H 
Pl*>([l'l']) = l, where 



P'»>([11]) = 
pl*>([l'l']) = 



1 



2(1 — cos 2 a) 
1 

2(1 — cos 2 a) 



3 Results 



The stability of a four-quark state can be ana- 
lyzed in terms of A E , the energy difference between 
its mass and that of the lowest two-meson threshold, 



A E =E 4q -E(M 1 ,M 2 



(7) 



where E Aq stands for the four-quark energy and 
E(M 1 ,M 2 ) for the energy of the two-meson threshold. 
Thus, A E < indicates all fall- apart decays are for- 
bidden, and therefore one has a proper bound state. 
A E > will indicate that the four-quark solution cor- 
responds to an unbound threshold (two free mesons). 
Thus, an energy above the threshold would simply 
mean that the system is unbound within our varia- 
tional approximation, suggesting that the minimum 
of the Hamiltonian is at the two-meson threshold. 

Another helpful tool analyzing the structure of a 
four-quark state is the value of the root mean square 
radii: (x 2 ) 1/2 , (y 2 ) 1/2 , and (z 2 ) l/2 . Bound four-quark 
states can be distinguished from two free mesons by 



| PQ | ) + ( | QP | * 



*|Pg|*)+(*|QP|* 



means of their root mean square radius 



RMS 



4g(2 9 ) 



4(2) ^<(, 



-R) 2 } 



1/2 



and in particular, their corresponding ratio, 



RMS 



RMS 



-RMS 



(6) 



(8) 



(9) 



iU 2 



where RMS M i+RMS A f 2 stands for the sum of the radii 
of the mesons corresponding to the lowest threshold. 

We show in Table[JJthe results obtained for several 
different four-quark states in the bottom and charm 
sectors. One can see how independently of their bind- 
ing energy, all of them present a sizable octet-octet 
component when the wave function is expressed in 
the (|lb|) coupling. Let us first of all concentrate on 
the two unbound states, A E > 0, one with S T = and 
one with S T = 1, given in Table [TJ The octet-octet 
component of basis (lb) can be expanded in terms of 
the vectors of basis (lc) as explained in the previous 
section. Thus, once expressions (6) are considered 
one finds that the probabilities are concentrated into 
a single physical channel, MM or MM*. In other 
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words, the octet-octet component of the basis (lb) 
or (lc) is a consequence of having identical quarks 
and antiquarks. Thus, four-quark unbound states are 
represented by two isolated mesons. This conclusion 
is strengthened when studying the root mean square 
radii, leading to a picture where the two quarks and 



the two antiquarks are far away, (x 2 ) 1 ^ 2 3> 1 fm and 
{y 2 ) 1 ^ 2 3> 1 fm, whereas the quark-antiquark pairs are 
located at a typical distance for a meson, (z 2 ) 1 ^ 2 < 1 
fm. 



Table 1. Four-quark state properties for selected quantum numbers. All states have positive parity and total 
orbital angular momentum L = 0. Energies are given in MeV and distances in fm. The notation M1M2 \t 
stands for mesons Mi and M2 with a relative orbital angular momentum £. P[|33) J 2 (|66) J 2 )] stands for the 
probability of the 33(66) components given in Eq. (|la[) and P[ ll) c ( | 88) c )] for the 11(88) components given 
in Eq. (|lbfl. Pmm, Pmm* , and Pm* m* stand for the probability of two pseudoscalar, pseudoscalar- vector or 

t.wn vprtnr tnpsntis 



(S T , I) 
Flavor 



(0,1) 



(Id) 



(1,0) 



(1,0) 
bbnn 



(0,0) 
bbnn 



Energy 
Threshold 



3877 
+5 



3952 
DD* | s 
+15 



3861 

DD* |s 
-76 



10395 
BB*\ S 
-217 



10948 
B\B \ P 
-153 



P[|33>«] 
P[|66>i 2 ] 



0.333 
0.667 



0.333 
0.667 



0.881 
0.119 



0.974 
0.026 



0.981 
0.019 



P[\n) c 

P[|88) c 



0.556 
0.444 



0.556 
0.444 



0.374 
0.626 



0.342 
0.658 



0.340 
0.660 



Pmm 
Pmm* 
Pm*m* 



1.000 



0.000 



1.000 
0.000 



0.505 
0.495 



0.531 
0.469 



0.254 



0.746 



(3.2)1/2 

<J/ 2 ) 1/2 
{z 2 ) 1 / 2 

RMSiq 



60.988 
60.988 
0.433 
30.492 
69.300 



13.804 
13.687 
0.617 
6.856 
11.640 



0.787 
0.590 
0.515 
0.363 
0.799 



0.684 
0.336 
0.503 
0.217 
0.700 



0.740 
0.542 
0.763 
0.330 
0.885 



Let us now turn to the bound states shown in Ta- 
ble [JJ A E < 0, one in the charm sector and two in 
the bottom one. In contrast to the results obtained 
for unbound states, when the octet-octet component 
of basis (lb) is expanded in terms of the vectors of 
basis (lc), equations (6) indicate a picture where the 
probabilities in all allowed physical channels are rel- 
evant. It is clear that the bound state must be gen- 
erated by an interaction that it is not present in the 
asymptotic channel, sequestering probability from a 
single singlet-singlet color vector from the interaction 
between color octets. Such systems are clear exam- 
ples of compact four-quark states, in other words, 
they cannot be expressed in terms of a single physical 
channel. Moreover, as can be seen in Table [TJ their 
typical sizes point to compact objects 20% smaller 
than a standard two-meson system. 

We have studied the dependence of the probabil- 
ity of a physical channel on the binding energy. For 
this purpose we have considered the simplest system 
from the numerical point of view, the (S T ,I) = (0,1) 
cenn state. Unfortunately, this state is unbound 



for any reasonable set of parameters. Therefore, we 
bind it by multiplying the interaction between the 
light quarks by a fudge factor. Such a modification 
does not affect the two-meson threshold while it de- 
creases the mass of the four-quark state. The re- 
sults are illustrated in Fig. [TJ showing how in the 
A E — > limit, the four-quark wave function is almost 
a pure single physical channel. Close to this limit 
one would find what could be defined as molecular 
states. Moreover, the size of the four-quark state 
increases when A E — > 0. When the probability con- 
centrates into a single physical channel (Pmm 1) 
the size of the system gets larger than the sum of two 
isolated mesons ^. We have identified the subsys- 
tems responsible for increasing the size of the four- 
quark state. Quark-quark ((cc 2 ) 1 / 2 ) and antiquark- 
antiquark ((y 2 ) 1 ^ 2 ) distances grow rapidly while the 
quark-antiquark distance {(z 2 ) 1 ^ 2 ) remains almost 
constant. This reinforces our previous result, point- 
ing to the appearance of two-meson-like structures 
whenever the binding energy goes to zero. 
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Fig. 1. Pmm as a function of A E . 

Although the present analysis has been performed 
by means of a particular quark interacting poten- 
tial ^, the conclusions derived are independent of 
the quark-quark interaction used. They mainly rely 
on using the same hamiltonian to describe tensors 
of different order, two and four-quark components 
in the present case. When dealing with a complete 
basis, any four-quark deeply bound state has to be 
compact. Only slightly bound systems could be con- 
sidered as molecular. Unbound states correspond to 
a two-meson system. A similar situation would be 
found in the two baryon system, the deuteron could 
be considered as a molecular-like state with a small 
percentage of its wave function on the AA chan- 
nel, whereas the H— dibaryon would be a compact 
six-quark state. When working with central forces, 
the only way of getting a bound system is to have 
a strong interaction between the constituents that 
are far apart in the asymptotic limit (quarks or an- 
tiquarks in the present case). In this case the short- 
range interaction will capture part of the probability 
of a two-meson threshold to form a bound state. This 
can be reinterpreted as an infinite sum over physical 
states. This is why the analysis performed here is so 



important before any conclusion can be made con- 
cerning the existence of compact four-quark states 
beyond simple molecular structures. 

If the prescription of using the same hamiltonian 
to describe all tensors in the Fock space is relaxed, 
new scenarios may appear. Among them, the inclu- 
sion of many-body forces is particularly relevant. In 
Ref. ™ the stability of QQnn and QQnn systems was 
analyzed in a simple string model considering only a 
multiquark confining interaction given by the mini- 
mum of a flip-flop or a butterfly potential in an at- 
tempt to discern whether confining interactions not 
factorizable as two-body potentials would influence 
the stability of four-quark states. The ground state 
of systems made of two quarks and two antiquarks 
of equal masses was found to be below the dissoci- 
ation threshold. While for the cryptoexotic QQnn 
the binding decreases when increasing the mass ratio 
mQ/m n , for the flavor exotic QQnn the effect of mass 
symmetry breaking is opposite. Others scenarios 
may emerge if different many-body forces, like many- 
body color interactions ^ or 't Hooft instanton-based 
three-body interactions * , are considered. 

4 Conclusions 

We have discussed the formalism to express the 
wave function of a general four-quark state in terms 
of physical channels, i.e., those constructed by using 
color singlet-singlet states. We have studied charmo- 
nium {ccnn) 1 bottomonium (bbnn) and exotic (ccnn 
and bbnn) four-quark states by means of a standard 
non-relativistic quark potential model. We look for 
possible bound states. Among them we are able 
to distinguish between meson-meson molecules and 
compact four-quark states. The importance of per- 
forming a complete analysis of the system, energy and 
wave function, in the vicinity of a two-meson thresh- 
old has been emphasized. 
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